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Consider a medium exhibiting memory. Such a situation arises in describing the structure of a solid far from 

■^ • thermodynamic equilibrium: in amorphous materials, under structural relaxation of high-T c oxide superconductors, 

under plastic deformation and destruction of solids, in solid solutions and in martensite macrostructures, and so on 

|[j. The existence of memory implies that if a force fit') acts on a system at an instant t' > 0, then a flux J arises 

O ' whose magnitude at a subsequent instant t > t' is given by 
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J(t) = I M{t - t')f(t')dt'. (1) 

For a system without memory, the time dependence of the memory function M(t — t') has the form 

M(t - t') = 7 5(t - t'), (2) 

where 7 is a positive constant and 8(t — t') is the Dirac (5-function. On substituting (2) into (1), we obtain the relation 

g- J r (t)=7/(*). (3) 

according to which in the absence of memory only a force value at the same instant t has an effect on the flux J(t). 
With memory inclusion, the 5-function in (2) spreads into a bell-shaped function with a width determining an interval 
r during which the force / has an effect on the flux J. For systems with ideal memory, we have r — > 00, i.e., the flux 
J(t) is formed over all the course of the action of the force /(*') up to the instant t. Formally, this is expressed by 
specifying the kernel of the integral relation ( 1 ) in the form 
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M(t - = 7 /t. (4) 

Here, the dependence on the instant t' of the force / action is absent, and the dependence of the memory on the flux 
f^v J measurement time t is given so as to satisfy the normalization condition 

f M(t - t')dt' = 7. (5) 
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Relation (1) as written in the time domain, is inconvenient because of the convolution (integral over t'). This can 
be eliminated by using the Laplace transformation 
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f{t) = J Z( A ) e ~ At ^-' /W = J o /(*)e Ai dt, (6) 

which enables one to move from the time t to the complex frequency A. After transforming both sides of the definition 
(1), this converts to the algebraic form 

J(A) = M(A)/(A). (7) 

Obviously the Laplace transform of the kernel (2), which corresponds to the absence of memory, reduces to a constant: 

M (A) = 7. (8) 

For ideal memory, we obtain from (4), in the limit |A|i 3> 1 
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M (A) = 7 /Ai. (9) 

Thus, the inclusion of memory leads to the transformation of the constant kernel (8) into the hyperbolic form (9). 

In this connection, one could ask: How does the function M (A) stand up if memory is complete but nonideal? This 
means that the memory manifests itself within the interval (0, t) preceding t but not at all instants t' . Assume, for 
example, that memory holds only at the points of a Cantor set. One can then expect that its fractal dimension D 
will be associated with a measure of memory preservation. 

Although elementary for simple time dependences (2) and (4), the problem of finding the Laplace transform for the 
memory acting at the points of the Cantor set is much more difficult. The appropriate calculations J2J arrive at the 
result 

M(A) = A D z- D , z=(l-£)Ai. (10) 

Here, the fractal dimension D is defined by usual equality p) D = In j/ln£ _1 , j is the branching index of the corre- 
sponding hierarchical tree, £ is the similarity parameter, and the constant Ajj = 2~ D ' 2 for the Cantor discontinuum 
with the number of elementary blocks j '• = 2. 

It can easily be seen that the memory function (10) satisfies the similarity condition M (£A) = £~ D M(A) with the 
exponent D. Being the fractal dimension of the Cantor set at whose points memory is switched on, this exponent 
thus represents a quantitative measure of the manifestation of memory effects. For an empty Cantor set (£ = 0), one 
has D = 0, and the dependence (10) reduces to the constant (8), corresponding to the entire absence of memory. The 
exponent D increases with the similarity parameter £ > 0, and the Laplace transform (10) becomes an increasingly 
more rapidly varying function. The limiting value £ = j -1 of the similarity parameter yields the dimension D = 1, 
which corresponds to the ideal memory governed by the function (9). 

Thus, systems with residual memory are described by the Laplace transform (10), where the exponent value 
< D < 1 determines the extent of memory preservation. Using the inverse Laplace transformation (6), we find for 
the time-dependent memory function in the simplest case of the Cantor set with j = 2 (see [|J , || ) 

M(t - t>) = (7/tt) [V2(l - £t] ~° T(l -D)(t- t') D -\ (11) 

where T(l — D) is the gamma- function. Substituting this into equality (1), we bring this equality in the form 

J(t) = (7/tt) [V2(l - £)] ° D- D f(t), (12) 

where the integral is introduced 

D- D f(t)=T(l-D) f (l-u^'/MdB 
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= D- 1 r(l-£>) f f((l-u)t)du D , 
Jo 



(13) 



whose fractional nature is reflected by the presence of the exponent D in the argument u of differential H . 

In the foregoing discussion we have used an integral representation for memory effects. It is easy to show, however, 
that an equivalent differential equation can be written that corresponds to this integral representation. For example, 
if we have some conserved quantity n (like the number density of atoms), then it space-time dependence n(r, t) is 
governed by the continuity equation 

m - - VJ < (14) 

where the flux J is given by equality (1). In the absence of memory effects, equality (3) is satisfied, where the force f 
is given by 

f = -V/x, (15) 

where /i is the chemical potential, and (14) leads to a conventional equation of diffusion type [|| 

^ = V( 7 V M ). (16) 



According to the earlier discussion, when memory is switched on, a fractional integral appears on the right-hand 
side of (16). The treatment of the resulting integro-differential equation, with partial derivatives and furthermore, of 
fractional order, presents a very difficult problem. 

However, this problem can be simplified when it is considered that the memory function M(t — t') vanishes within 
some intervals, thus leading to a fixing of the flux J in (14) and correspondingly to the vanishing of the velocity dn/dt. 
Therefore, it would appear natural that loss of memory can be taken into account by a fractional property, not only 
of the flux integral (1) but also of the time derivative in (14). Because equation (16) corresponds to the entire absence 
of memory (D = 0), when it is switched on, as is reflected by the growth of the exponent D, the rate of n variation 
must decrease, and it is reasonable to postulate a fractional-differential equation of the form (see further) 

d q n 

— =V( 7 Vm), q=l-D. (17) 

As would be expected, the growth of the memory parameter D leads to lowering the order q of this equation. In 
systems with ideal memory (D — 1), the derivative on the left-hand side of (17) disappears so that one should treat 
the behavior of the flux (1) itself which is what we have done above. 

Our considerations in favor of the specific form q = 1 — D for the fractional degree in equation (17) act as guidelines. 
To confirm these, we will show that equation (12) for the flux in the form of the fractional integral (13) is in agreement 
with the solution of the fractional-order differential equation (17). Indeed, replacing the variable u by (1 — u)t in 
the integrand of the second integral of (13), one can easily show the time dependence of the flux to have the form 
J ex t~ D . Substituting this dependence into (14), we find the solution 

noct 1 ^, (18) 

satisfying the tested equation (17). This confirms conclusively the equivalence of the mutually complementary concepts 
of the fractional integral (13) and the fractional-differential equation (17). Notably is displayed in the same form the 
geometric relation L oc l 1 ^ for coastline length L{1) and in the physical dependence n(t) obtained. In particular, it 
follows that the time t plays the some role as the length of the segments I, covering a D-dimcnsional fractal set in the 
space of states, the points of this set determining the memory of the system. 

When interpreting equation (17) in fractional derivatives, we start from the fact that in the usual case D = this 
equation governs irreversible processes like diffusion, heat conductivity, and so forth y. In the latter processes, the 
microscopic memory with respect to time inversion is completely lost: the mechanical equations of motion for an 
isolated object exhibit ideal memory, manifesting itself in the invariance with respect to the replacement t — > — t, 
whereas such invariance is completely violated in the thermodynamic equation (16). It follows that a decrease in the 
derivative order q = 1 — D in (17) corresponds to the switching on of memory channels whose part is determined by 
the fractal dimension D. The remaining q = 1 — D channels provide the irreversibility of the system. In the initial 
form, the residual irreversibility can be taken into account by using the fractional form of the Liouville equation [|| 

ih d D p rTT . , „. 

-— % = [H,p], 19 

t au u 

where Ti is the Planck constant, H is the Hamiltonian, p is the nonequilibrium statistical operator, t is the evolution 
time of a macroscopic system, and u = t'/t is the dimensionless microscopic time bounded by the condition u < 1. 
It is apparent that the condition D < 1 enables one, in the framework of (19), to take into account the irreversible 
effects in the interval (0, t) which are due to the loss of q = 1 — D deterministic channels. Such a means of allowing 
for the irreversibility is obviously preferable to the phenomenological addition to the right-hand side of (19) of the 
relaxation term —p/r, where r is the relaxation time, or to the inclusion of an infinitesimal source B). 

Apart from equations of the diffusion type, one might expect the appearance of fractional derivatives in a case of 
the description of the motion of particles under inelastic collisions. Performing calculations similar to those that led 
to the flux of the form (12), one can show [g] that, if a force F acts on a particle of mass m in each collision, then the 
change in its velocity 
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2(1 - [t/m)D- D F (20) 



is determined by a fractional integral of the form (13). From the above example of the diffusion equation, one can 
see that, in order to switch to the corresponding fractional-differential equation, it is sufficient to operate on equation 
(20) with the operator D D = d D jdvP which is inverse with respect to the fractional integral D~ D . In the case of an 
elastic force F = — Af2V 2 r, where A is the elastic constant and Vt is the atomic volume, we then obtain the generalized 
transport equation for the particle coordinate r H 



d 1+D r 

—jy + CrfjV^O, (21) 

where the dimensionless time u = t 1 /t < 1 is used and the characteristic velocity c, defined by 



c 2 = (Aft/m) [\/2(l - £) 
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(22) 



is introduced. Equation (21) governs a new type of wave motion that is intermediate between conventional diffusion 
(D = 0) and classical waves (D = 1). Accordingly, formula (22) determines the diffusion coefficient D — c 2 t in the 
first case and the wave velocity c in the second case. 

The fractional-differential equations (17) and (21) are written as applied to the study of the space-time behavior 
of conserved quantities. As we known ||, this is reflected in the presence of a second derivative with respect to the 
coordinate in the right-hand sides of these equations. For a nonconserved quantity 77, this derivative vanishes, and 
the diffusion equation (17) converts to the generalized Landau-Khalatnikov relaxation equation 

Jb=S = -to*)* (23) 

and the wave equation (21) converts to the generalized equation for an harmonic oscillator 
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V + (uotf^ = 0; (24) 



here, as above, u = t'/t is the dimensionless time, and the parameters 725, ojd oc [v2(1 — £)] determine the 
relaxation time to = 7^ and the oscillator frequency u>d- The fractal dimension D specifies a measure of the 
residual memory of the system. In its absence (D — 0), both equations (23) and (24) reduce to the equation of 
Debye's relaxation which governs the behavior of the simplest thermodynamic systems (70 = to^t in this case). For 
ideal memory (D = 1), equation (23) degenerates into the condition of phase equilibrium rj = r] and (24) degenerates 
into the ordinary oscillator relation. 

The above consideration shows that, in limiting cases, one can recognize the following equations of motion, which 
correspond to the different types of behavior of a continuous medium. First, there is the continuity equation, being the 
condition for the conservation of particles in the medium. In the context of a mechanistic approach, the oscillatory 
behavior is governed either by an equation of oscillation, which involves a second derivative with respect to time, 
or by a wave equation involving second derivatives with respect to both time and space coordinates. Under the 
thermodynamic description of a system determined by a hydrodynamic mode, its amplitude (order parameter) may 
be characterized both by reactive and by dissipative regimes. In the first regime, the equation of motion reduces to 
an oscillatory equation for a nonconserved parameter and to a wave equation for a conserved one. Upon entering the 
dissipative regime, the time-derivative order is lowered to the first order. 

Although these equations govern quite different physical situations, our discussion shows them to transform smoothly 
from one to another H, W. Thus, the question arises: Is it possible to give a model of the medium whose equation 
of motion includes all the above types as special cases? Following Iq], we will show how such an equation can be 
obtained. 

We start from the continuity equation 77+ Vj = for the order parameter field r/(r, t) characterized by the flux j(r, t). 
In the general case, the latter is nonlocally related to the distribution of the chemical potential /u(r, t). In the framework 
of the Fourier representation with respect to the space coordinate and of the Laplace representation with respect to 
time (we take into account the causality condition), this relation reduces to the equality j(k, A) = M(k, A)/i(k, A), 
where k is the wave vector and A is the complex frequency. 

To find the dependence of the memory function M(k, A) on k, we will consider the simplest one-dimensional model 
of the medium in the form of a segment < x < I. Then the flux through the right boundary is given by (cf.Eq.(l)) 

j(l) = / M(l - x)n(x)dx. (25) 



We assume that the function M(x) takes on nonzero values only at the points of the Cantor set which is obtained 
from the original segment [0,^] by the convolution method. Then, as in Ref. |2j], it can be shown that the desired 
dependence is given by the equality similar to (10): 



M(k)=&- aU KT u , K = ikl(l-0- (26) 



Here, the constant a ~ 1 is introduced (according to 0, a = 2~ 1 ln2 = 0.347 for the simplest Cantor set), D, £, and 
I are the fractal dimension, the similarity parameter, and the scale, respectively, which describe the set of points on 
the 2-axis where the order parameters conserved. Using (26), we find the dispersion law from the continuity equation: 

\ = B_ D k 1 - D 1 B^ D =i 1 - D c- aD {l-S,)- D l l - D d t i/di 1 . (27a) 

When we reduce the fractal set to a continuum (£ = 1/2), its dimension is then D = 1, and relation (27a) corresponds 
to the Landau-Khalatnikov equation for the nonconserved order parameter. In the general case of £ < 1/2, upon 
reduction of the fractal dimension of the set with ordering points, the vanishing of long- wave harmonics occurs, which 
is inherent in the order parameter conserved. However, even the maximum value of the exponent 1 — D, corresponding 
to the empty set (£ = 0), proves to be half as large as the exponent 2 inherent in the order parameter conserved. 

To raise the exponent in the dispersion law (27a), we assume that the order parameter becomes conserved after going 
to the empty set as a result of the reducing the similarity parameter £ from 1/2 to 0. Next, using a procedure inverse 
to convolution, we increase the fractal dimension to a finite value D > 0. Obviously, this can be done by decreasing 
the similarity parameter £ from oo to a finite value £ > 2. The condition £ > 1 implies that we stretch the segments 
rather than convolute them. Then, the functional equation M(k/£) = 2~ 1 M(k) takes the form M(£k) = 2M(k). Its 
solution M(k) — e aD k d has a positive exponent whose value is now determined by the relation 

L» = ln2/ln£ (28) 

instead of the usual equation D = — ln2/ln£. Accordingly, the dispersion law for the partially conserved order 
parameter has the form 

X = B D k l+D , B D =i 1+D e aD \l-£\ D l 1+D dn/dr]. (27) 

As seen from (27a), the distinction from the case of the conserved order parameter reduces to the change of signs 
of D and £ — 1. According to (28), this corresponds to the transition from contraction (£ < 1) to stretching (£ > 1) 
when constructing a fractal. 

Thus, if it is assumed that the similarity parameter £ is defined in the regions (0, 1/2] and [2, oo), the dispersion law 
will encompass both the cases of the nonconserved (D < 0) and conserved (D > 0) order parameter. The transition 
between the limiting cases occurs due to a smooth variation in the similarity parameter £. At £ = 1/2, there is 
D = — 1, and we have complete nonconservation of the order parameter. A decrease in £ down to values < £ < 1/2 
partially violates the nonconservation conditions, thus leading to an increase in the fractal dimension in the range 
— 1 < D < and to the rise in the exponent of the dispersion law (27). With the empty Cantor set (D = 0), we 
go from the nonconserved order parameter (£ — v 0) to the conserved one (£ — > oo). A decrease in the similarity 
parameter within the segment 2<£<oo(0<Z?<l) provides the rise in the number of conservation channels, thus 
increasing the exponent of the dispersion law (27). Finally, with the complete conservation of the order parameter 
(£ = 2, D = 1), we obtain, as expected, A oc k 2 . 

It is easy to understand that this variation in the order of the derivative with respect to the space coordinate 
can be extended to the derivative with respect to time. Indeed, as shown in pi, the memory switching on relative 
to the time inversion t — > — t in the equation of particle motion, leads to the appearance in (27) of the exponent 
1 — Dq at frequency A, where Dq is the dimension of a fractal set on the time axis, at the points t of which the 
mechanical memory is switched on. The set itself is represented by the similarity transformation F(z/£) = 2~ l F(z), 
where z — Ai(l — £o) is the argument and < £o <• 1/2- The transition from the contraction region £ <• 1/2 into 
the stretching region £o > 2, corresponding to the functional equation M(t;oz) — 2M(z), occurs on the empty set 
(Dq = 0) and implies the replacement of the corpuscular memory by the wave memory The latter is enhanced with 
a decrease in the similarity parameter £o and becomes ideal at £o = 2. Then, the exponent at the frequency A in (27) 
takes on the value 2 inherent in the wave equation. 

Thus, in the general case of partial nonconservation/conservation of the order parameter and of the non-ideal 
memory of particle/wave, the dispersion law for a one-dimensional fractal medium takes the form 

q=(l + D)/(l + D ), 

c aD \£-l\ D l 1+D d[i 



eaoA>|£ -l|Aj r i-kDo dry' 



ao,a~l, (29) 



where / and r are the characteristic space-time scales, and the dependence of the dimension D upon the parameter 
£o is determined by the same equality (28) as for the space component £>(£)• Accordingly, the fractional-differential 



equations of motion are written in the forrnn 

dt l + D + % dx l + D - U - ^) 

In the linear approximation /i oc 77, the constant c may be factored out of the differentiation sign. 

It can be easily seen that, in the limiting cases of the integer-valued space-time dimensions D Q and D, the linearized 
equation (30) encompasses the main regimes of medium evolution. For example, in the case Dq — 1, this governs the 
wave motion at D = I , uniform oscillations at D = — 1, and localized oscillations at D — 0. The constant c equals the 
square of wave velocity, the square of oscillation frequency, and the product of localization region size and oscillation 
frequency, respectively. In the case Dq = 0, we have the Landau-Khalatnikov equation for the conserved (D = 1) 
and nonconserved (D = — 1) order parameters and the traveling- wave equation (D = 0). The constant c reduces to 
the diffusion coefficient in the first case, to the inverse relaxation time in the second case, and to the wave velocity 
in the third case. And finally, in the static case Do = — 1, we obtain the Poisson equation at D = 1, the equilibrium 
condition at D = —I, and the Debye screening equation at D = 0. Here, c reduces, respectively, to the square of 
screening radius, to the inverse susceptibility d^/dr], and to the screening radius. 

It is evident that the approach outlined enables one to represent not only linear equations of motion but also to 
generalize these to the nonlinear case and to take into account the higher dispersion terms. For example, if, by analogy 
with the theory of solitons M, we set that, in addition to the term ex k D , the flux involves the term ex k 2+A with the 
exponent — 1 < A < 1, then a dispersion component of a higher order 3 + A appears. The nonlinearity is reflected 
in the dependence fJ.(j)) of the chemical potential on the order parameter. One should take into account that the 
dimension S of a nonlinear channel does not reduce to the corresponding quantity D for a linear one. As a result, 
equation (30) takes the form 

d 1+I S , . 1+D d^ D V d 3 +% d^/fo) _ 

d^+Do C °d.T 1 + D ' da; 3 + A dx^ 6 



e^fo-lfl^ 5 dj 



m = „, n ,r „ n „ , n :£: (3i) 



where Co corresponds to rj ~ 0, > is the dispersion parameter, and ji(r)) is the nonlinear contribution into the 
chemical potential. This equation reduces to the Korteweg-de Vries equation at Dq = D = A = S = and / ex i] 2 , 
to the nonlinear Schrodinger one at Do = D = 0, A = S = — I, / = rf, and finally to the sine-Gordon equation at 
Do = D = 1, /? = 0, S = — 1,/ = sin 7/. As is known, these equations are notable in that their symmetry within 
the framework of the inverse problem of the scattering theory provides the existence of soliton solutions |9| . Here the 
question arises: Will this symmetry also manifest itself at other values of the differentiation orders D, Do, A, and 
6 and/or with different choices of (3 and f(j])7 To answer this question demands investigation of the symmetry of 
multidimensional fractal sets and, to our knowledge, this has not been done. 

Now, one can note only that conformal symmetry (see Jl0|] ) states the follow connection between the above expo- 
nents: 

D = D Q = 2(1 + 5 -n) -1 . (32) 

Here exponent n determines the nonliniar contribution f(r/) oc rf 1 ^ 1 . Within framework of the Landau ?7 4 -theory, one 
has n — 4 for a thermodynamic system characterised by the additive noise which intensity is the temperature. But 
with passing to the multiplicative noise with the amplitude proportional to 77°, < a < 1 the powers of rj for all terms 
in Eq.(3f), excluding the first, are lowered by the value 2a |ff|. 



1 In the d- dimensional case, the second term of equation (30) is replaced by a sum of terms characterized by scale h, fractional 
dimension Di, similarity parameter f, and constant a;, for i — 1, . . . ,d. 
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